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Of NBOs and kHz QPOs: a low-frequency modulation 
in resonant oscillations of relativistic accretion disks 
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Abstract 

The origin of quasi periodic modulations of flux in the kilohertz range (kHz QPOs). observed in low- 
mass X-ray binaries, is usually assumed to be physically distinct from that of the "normal branch oscilla- 
tions" (NBOs) in the Z-sourccs. We show that a low-frequency modulation of the kHz QPOs is a natural 
consequence of the non-linear relativistic resonance suggested previously to explain the properties of the 
high-frequency twin peaks. The theoretical results discussed here are reminiscent of the 6 Hz variations of 
frequency and amplitude of the kHz QPOs reported by Yu, van der Klis and Jonker (2001). 
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1. Introduction 

Power density spectra of the X-ray flux of low-mass 
X-ray binaries have a rich phenomenology almost en- 
tirely lacking a physical explanation (van der Klis 1989, 
2000). Twin kHz QPOs in neutron star sources, and 
their kHz counterpart in black-hole systems, have at- 
tracted the most attention because their frequency is com- 
parable to the orbital frequencies in the innermost dy- 
namically stable regions of the accretion flow (Kluzniak, 
Michelson & Wagoner 1990; van der Klis et al. 1996). 
Ad hoc models of QPOs have not withstood the test 
of time, but a more fundamental approach of study- 
ing the eigenmodes of accretion disks remains promising 
(Wagoner 1999; Kato 2001). It seems necessary to include 
non-linear effects to be in agreement with observations 
(Kluzniak & Abramowicz 2003). Here, we wish to point 
out that even the simplest non-linear equations of motion 
have a rich structure which is yet to be explored fully in 
the context of QPOs. 

It is now recognized that non-linear coupling in the 
motion of accreting fluid may be responsible for some of 
the observed features of (kHz) QPOs in neutron-star and 
black-hole systems. It has been suggested that a non- 
linear resonance may be responsible for these twin QPOs 
(Kluzniak & Abramowicz 2001; Kato 2003; Kluzniak & 
Abramowicz 2003). Among the reasons to believe that 
these highest frequencies reflect a non-linear resonance be- 
tween two oscillation modes of a disk, probably occurring 
only in strong-field gravity, are the 3 : 2 ratio of the fre- 
quencies of the twin QPOs in black holes (Abramowicz & 
Kluzniak 2001; McClintock & Remillard et al. 2004), and 
the sub-harmonic frequency difference between the two 
QPOs observed in the accreting 2.5 ms pulsar (Wijnands 
et al. 2003; Kluzniak et al. 2004; Lee, Abramowicz & 



Kluzniak 2004). It has been found that a parametric-like 
resonance can be excited between the radial and vertical 
epicyclic frequencies in a simplified mathematical model of 
fluid motion in Einstein's gravity, and that the ratio of the 
two kHz frequencies can be made to match closely the one 
observed in Sco X-1 if the source is slightly off-resonance 
(Abramowicz et al. 2003; Rebusco 2004). In this paper, we 
employ the same model of twin kHz QPOs and we inves- 
tigate its consequences for the time behavior of predicted 
frequencies and amplitudes. We lay down our description 
within the framework of an (internal) resonance model 
describing the oscillatory modes of a conservative system. 
Details of the adopted method and calculations can be 
found elsewhere (Horak 2004). 

2. The model 

We solved a very general system of two coupled nonlin- 
ear oscillators described by the following governing equa- 
tions near the resonance : tjj. ~ 3 : 2, 

5r + ujl6r = ujlf{6r,5e,6r,5e), (1) 

5e + ujl5e = L4g{5r,6eJrJe), (2) 

under the condition that these two equations are invariant 
with respect to reflection of time, i.e., the Taylor expan- 
sion of functions / and g docs not contain odd powers of 
the time derivatives of 5r and 59. These equations may be 
taken to approximate two modes of an accretion disk with 
eigenfrequencies close to the radial and vertical epicyclic 
frequencies, tHr and ug, respectively. The functions /, g 
depend on the space-time metric, as well as on the prop- 
erties of the fluid flow. Any particular mechanism should 
specify the physical meaning of these functions and fix 
their values, e.g. in terms of the accretion rate, resonance 
radius and other parameters. Such special models have 
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been indeed proposed and they are encompassed by the 
scheme developed herein. For example, a special form of 
/ and g was assumed in Abramowicz ct al. (2003) and 
Rebusco (2004), based on the idea of parametric reso- 
nance, but here we keep the discussion general and use 
the special solution for the purpose of an example and a 
numerical check, as described below. In fact, from the 
mathematical point of view, the free parameters of the 
model are the expansion coefficients through the fourth 
order of the functions / and g, which unambiguously de- 
fine the solution. 

Wc applied the method of multiple time scales (Nayfeh 
& Mook 1979) and looked for a solution in the form of a 
series. 



5r{t) = eri 
5e{t) = eOi 



0{e' 

e'02 + e^e3 + e^^4 + O(e^ 



(3) 
(4) 



This approach enables us to find a uniformly converg- 
ing solution, which does not suffer from secularly growing 
terms. We introduced new independent variables Tk = e^t 
and derived the conditions for the elimination of secular 
terms in the "fast" variable Tq = t, up to the fourth order 
in the expansion parameter e ^ 1. We restrict ourselves to 
the leading terms of the expansions, which take the form 



5r{t)^Ar{t)e''^-'- 
6e{t)^Ag{t)e"^''' 



' c.c. 



(5) 
(6) 



The solvability conditions describe the time evolution of 
complex amplitudes, which wc further rewrite in the form 



J4>r 
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Ae = -age 



(7) 



The amplitudes and phases therefore satisfy equations 



■a^al sin7. 



P^e 3 

-— ae sm7, 
Id 



-— [Kra^ + Keag\ 



Xean 



■ ar ag C0S7, 



-a;; C0S7, 



(8) 

(9) 
(10) 

(11) 



aujr 

Te 

f3uje_ 
16 

where 7 = 2(j)e — 3(f)r — crt and a = Scj,, — 2uje are a phase 
function and a detuning parameter, and a, /3, Kr, ne, Ar 
and Xg are real constants describing our system. The 
detuning parameter a describes small departures of the 
eigenfrequencies from sharp rational ratio. In this way, 
the method of multiple scales allows us to capture the 
amplitude-frequency interaction that is present in the os- 
cillating solution. 

Equations (10) and (11) provide frequency corrections 
Aujr = (f'r and Awg = (j)g to the eigenfrequencies Wr and 
ujg. The actual, corrected frequencies will be denoted by 
an asterisk in order to distinguish them from the corre- 
sponding eigenfrequencies: v* = {lu r + Aojr)/ {2Tr), Vg = 
{ujg + Au!g) / {2tt) . Mathematical validity of the adopted 
approximation requires the frequency corrections to be 
small in our model. 
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Fig. 1. Time evolution of amplitudes (top panel) and the 
epicyclic frequencies (vertical in the middle, radial at bottom) 
of high-frequency QPOs. The two epicyclic modes are in- 
terpreted as radial/ vertical oscillations of the accreting fluid. 
The variation of the oscillation amplitudes is not arbitrary in 
the adopted resonance model: instead, a relation between the 
amplitudes and frequencies is naturally predicted. 

Equations (8) and (9) can be solved for sin7 and inte- 
grated to read 



PuJg 



E = const. 



(12) 



The quantity E is proportional to the energy of oscil- 
lations and is conserved in this order of approximation. 
Equation (12) describes an ellipse in {ar,ag) coordinates 
and it rules the flow of energy between the oscillation 
modes. 

The two epicyclic frequencies ujr and ujg define the ba- 
sic modes of oscillations. Secular terms and the possibility 
of parametric resonance appear in the second order of ap- 
proximation. In the third order, the frequency corrections 
are characterized by the foursome of parameters: Kr, Kg, 
Xr and Xg. Finally, the two additional parameters, a and 
P, provide the fourth-order corrections, and describe the 
semi-axes ratio of the energy ellipse. The above-given list 
completes the most general description of any system de- 
scribed by equations of the form (l)-(2). 

The general solution corresponds to the periodic ex- 
change of energy between the two oscillators. If the sys- 
tem is far from the steady state, the amplitudes and fre- 
quencies of oscillations fluctuate, maintaining the energy 
condition (12). The period of energy exchange is approx- 
imately given by 



16^ £;-3/2_ 

f3ujg 



(13) 
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On the other hand, particular steady-state solutions are 
possible in which the observed frequencies of the two 
modes remain strictly constant in time in 3 : 2 ratio, de- 
spite the fact that the eigcnfrcqucncics depart from it. 
Such solutions play a role of singular points of the system 
(8)-(ll). Close to steady state the approximation (13) 
fails and the period becomes greater (Horak 2004). 

In order to verify accuracy of our work, we checked that 
the numerical solution of a special form (Abramowicz et 
al. 2003) of the system of equations (l)-(2), in which 
the right-hand side corresponded to the case of nearly 
geodesic circular motion in the pseudo-Newtonian poten- 
tial of Paczynski and Wiita, does indeed closely follow the 
ellipse of eq. (12). Therefore we can be confident that the 
above-described analytical method gives credible results. 

It is evident from eqs. (8), (9), (12), and (13), that the 
amplitudes and frequencies of high-frequency QPOs are 
modulated at the frequency 1/3 = l/T, which is related to 
the amplitudes of the original high-frequency QPOs. We 
remind the reader that here we restrict ourselves to the 
case We : = 3 : 2, however, similar results can be obtained 
also for other eigenfrequency ratios. 

3. Low-frequency modulation 

In the general discussion here, our model has six free 
parameters, corresponding to (a combination of) the low- 
est order expansion coefficients in the Taylor series of the 
unspecified functions / and g in equations (1) and (2). 
Hence, we cannot predict unique behavior. However, it 
is remarkable that a modulation of frequencies and am- 
plitudes follows naturally from the assumptions made. 
Because the oscillators are non-linear, their frequency 
varies with amplitude. Because the two oscillators arc 
coupled in a system with constant energy, the amplitudes 
of the oscillators arc anticorrclatcd. Because the two os- 
cillators are in resonance, their frequencies are correlated. 

It is not lost on us that similar correlations — discussed 
by Yu et al. (2001) — hold also on long timescales. For 
example, van der Klis et al. (1997) show that as the kHz 
frequencies increase, the ratio of power in the upper to 
lower kHz QPOs decreases (up to a point). Our Figure 1 
demonstrates this to be the case also here. We note 
that in previous work, a particular relation was found to 
hold between the two frequencies of the system, in agree- 
ment with that observed in long-term variations of the 
twin kHz QPO frequencies in Sco X-1 (Abramowicz et 
al. 2003; Rcbusco 2004). Here, we consider variations on 
shorter timescales. 

As we mentioned above [see equation (13)], the 
timescale of modulation in our model is directly related to 
energy i.e., the weighted sum of amplitudes squared, 
which we adjusted to obtain a modulation at the NBO 
frequency of Sco X-1, 1/3 w 6 Hz. We identify the corrected 
frequencies of our model, i.e., v* and Vg, with the twin 
kHz QPO frequencies t'lower and i^upper, respectively. 

Figure 1 (top panel) exhibits the time variation of the 
"radial" and "vertical" amplitudes of oscillations. In the 
lower panels we show the correlation between the two kHz 



frequencies found in our solution. By the assumption of 
nearly 3 : 2 resonance, the two frequencies of actual oscil- 
lations satisfy relation f upper ~ LSi^f^^^j.. Notice a perfect 
correlation between the variation of the lower amplitude 
and the variation of the upper frequency v'^p-p^r- This we 
interpret as the same correlation that was found by Yu 
et al. (2001); compare their figure 2. The magnitude of 
frequency variation agrees to within a factor of 3 with the 
data: it is 20 Hz for I'uppcr = 1.1 kHz in the data, and 
about 7 parts in 1000 in our calculation. 

4. Discussion 

We have found that a non-linear resonance in a system 
simulating an accretion disk, invoked previously to explain 
the appearance of two frequencies in approximately 3 : 2 
ratio in black-hole and neutron-star X-ray data, results in 
a periodic time variation of the frequencies. 

The Fourier transform of the frequency-modulated and 
amplitude-modulated periodic signal would result in sev- 
eral sidebands, in practice leading to an increase in width 
of the (noisy) signal. If the observed radiation flux were 
modulated with the squared modulus of the amplitudes 
(5)-(6), its Fourier transform would exhibit a (weak) com- 
ponent at low frequency of the modulation apparent in 
figure 1, in addition to the high-frequency signal with its 
sidebands. In this exploratory work, we are not attempt- 
ing to model the full power spectrum of QPO sources and, 
as yet, we have made no attempt to translate the ampli- 
tudes of motion in the model into modulations of the X- 
ray flux. Strictly speaking, our toy-model gives a coherent 
signal rather than a QPO, and no details of the excitation 
or damping were modelled. 

It has been suggested that the high-frequency QPOs 
vary on a timescale of seconds in some sources, notably 
in Sco X-1 (Yu et al. 2001) and in the black hole candi- 
date XTE J1550-564 (Yu, van der Khs & Fender 2002). 
The low-frequency modulation occurs at « 6 Hz in both 
sources. In the above described calculation we were able 
to reproduce this variation for Sco X-1, including the anti- 
correlation between the amplitude of the lower peak and 
the frequency of the upper one. A similar approach is pos- 
sible also in case of XTE J1550-564, although the black 
hole candidates typically exhibit lower frequencies com- 
pared to those in neutron stars, and particularly to those 
in Sco X-1. This means that if the 6 Hz QPO seen by 
Yu, van der Klis & Fender 2002 corresponds to the mod- 
ulation discussed in our model, then the ratio fiowor/^'3 
has to be set differently (about 30 in XTE 1550-564 with 
'-'lower = 184 Hz). And this in turn implies E is different 
in both systems. 

If the correspondence of our results with the observed 
modulation of kHz QPO properties on the 6 Hz NBO 
timescale is not accidental, for the first time we would 
have a physical explanation for the presence of the rather 
low frequency QPO in what is otherwise a domain of rapid 
variability. The larger point is that non-linear resonance 
(most likely between modes of oscillation possible only 
in strong gravity) holds promise for explaining not only 
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the highest frequencies observed in accreting neutron stars 
and black holes, but also the mysterious phenomenology 
of low frequency features in the power density spectrum, 
without invoking additional mechanisms. 
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